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CHARACTERIZATION OF ACOUSTIC DISTURBANCES 
IN LINEARLY SHEARED FLOWS 

S. P. Koutsoyannis 

Joint Institute for Aeronautics and Acoustics 
Department of Aeronautics and Astronautics 
Stanford University, Stanford, California 94305 

The equation describing the plane wave propagation, the stability 
or the rectangular duct mode characteristics in a compressible inviscid 
linearly sheared parallel, but otherwise homogenous flow, is shown to 
be reducible to Whittaker's equation. The resulting solutions, which 
are real, viewed as functions of two variables, depend on a parameter 
and an argument the values of which have precise physical meanings de- 
pending on the problem. The exact solutions in terms of Whittaker func- 
tions are used to obtain a number of known results of plane wave propa- 
gation and stability in linearly sheared flows as limiting cases in 
which the speed of sound goes to infinity (incompressible limit) or the 
shear layer thickness, or wave number, goes to zero (vortex sheet limit). 
The usefulness of the exact solutions is then discussed in connection 
with the problems of plane wave propagation and the stability of a finite 
thickness shear layer with a linear velocity profile. With respect to the 
plane wave propagation it is shown that, unlike the compressible vortex 
sheet, the shear layer possesses no resonances and no Brewster angles, 
whereas with respect to the stability problem it is shown that again unlike 
the compressible vortex sheet, the thin shear layer is unstable to long 
wavelength disturbances for all Mach numbers. These results imply that 
the reflection and stability characteristics of a nonzero thickness but 
thin shear layer (i.e., the long wavelength characteristics) do not go 
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over smoothly into the results of the compressible vortex sheet as the 
wave number approaches zero except for a limited range of generally sub- 
sonic relative flow of the two parallel streams bounding the shear layer. 

1. Introduction 

Although problems pertaining to plane wave propagation, stability 
and rectangular mode in a compressible inviscid linearly sheared paral- 
lel flow have received considerable attention, with the exception of the 
work of Goldstein and Rice 11], all previous work has been concerned 
with either asymptotic and/or series solutions of the governing equation 
which is of course the same for all three above classes of problems. 
Typical examples of earlier work are that of Kuchemann [2] and Pridmore- 
Brown [3], and of later work, that of Graham and Graham [4] and Goldstein 
and Rice [1]. Kuchemann [2] obtained formal series solutions of the den- 
sity perturbation equation and also asymptotic solutions purportedly valid 
for large Mach numbers. His series solution, although it is given in a 
cumbersome and lengthy form, is correct and agrees with our compact form 
(see equations (8) and (9) below) , but his asymptotic solution is in error 
(see discussion following equation (17) below) . Pridmore-Brown [3] solved 
the pressure perturbation equation in the short wavelength approximation. 
His asymptotic solution which is in terms of Airy functions is limited in 
that the interesting region r) ■* 0 is excluded. Graham and Graham [4] ob- 
tained a series solution for the density perturbation equation apparently 
unaware of the earlier work of Kuchemann. Their series solution agrees 
with that of Kuchemann and with our series expressions mentioned above. 
Goldstein and Rice [1] were apparently the first to obtain solutions to 
the governing equation in terms of special functions. Using an unusual 
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double transformation (the essential part being a differential transfor- 
mation w„th a gaussian kernel) they were able to obtain from the original 
second order equation for the pressure perturbation an exact third order 
equation from which two independent solutions of the original second or- 
der equation were easily extracted in terms of the parabolic cylinder D 
functions. Unfortunately the solutions obtained by Goldstein and Rice 
(1J are not in terms of single parabolic cy Linder D functions, but combi- 
nations of D functions of different orders and in addition the solutions 

obtained were complex. It is not possible to reduce the solutions given 

* 

by Goldstein and Rice (11 to our solutions which are in terms of single 
Whittaker M or W functions, which is not surprising since our solutions 
are real, whereas those of Goldstein and Rice are complex. But if one 
does choose either specific linear combinations of the solutions of Gold- 
stein and Rice or the so-called even and odd solutions of Weber's equation 
instead of the parabolic cylinder functions D used in the above work, then 
it may be shown that our real solutions, in terms of single Whittaker 
functions and those of Goldstein and Rice are compatible. Recently in a 
study Jones [5] considered the stability of an everywhere subsonic (M < 1) 
shear layer with a linear velocity profile with the result that for 
0 < M < 1 there is a characteristic Strouhal number below which the layer 
is unstable. The problem investigated by Jones [5] is substantially dif- 
ferent from the problem treated here. In the first place he considers a 
source at a finite position from the layer whereas we consider plane waves 
emanating from -°° in the propagation problem and the usual formulation of 
the eigenvalue equation for parallel flow stability (see Betchov and 
Criminale [6]) without the necessity of invoking causality which Jones 
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does, Moreover, Jones' work pertains to the open-ended region 0 < M < 1 
whereas our3 covers the whole range M J> 0 including the incompressible 
limit M * 0 and the supersonic flow regime. 

2. Basic Differential Equation and Its Solution 

In a homogeneous inviscid compressible parallel shear flow having 
a linear velocity profile in the z-direction only, i.e., U ■ U(z) * bz, 
with b constant, it may be easily shown that, starting either from the 
linearized equations of motion or directly using the appropriate linear- 
ized form of the convective wave equation, the z-dependent part p(z) of 
the pressure perturbation p(r) is governed by the equation: 

P nn - | P n + <4T) 2 (n 2 - l) P “ 0 

where 

T1 * ~ - M , M = M(z) = , 4t = ~ K (2) 

K a D 

M is the local Mach number, a the (uniform) sound speed, and K and u> ac- 
quire the following meanings depending on the problem at hand: 

For propagation of a plane wave of wave vector k and frequency 0) 
impinging on the shear layer from a half-space (z < 0) of relative 
rest and at an angle 0 measured from the z-axis ~ < 0 < + jj, 

K = sin 0 . (2a) 
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For the stability of a free shear layer to assumed disturbances 
of the form p(r) « p{ 2 )e iot ^ x ’ ct ^ (a and c possibly complex), 


K * 2. and 0 ) » ac . (2b) 

C 

For sound propagation in rectangular ducts and for assumed distur- 
bances of the form p(r) « p(z)e ia <kx ~ Wt * (k a nd w real), 


K * — 

to 


(2c) 


The transformation 


1 

p «• n 2 w(?) , c ■ qn 2 O) 

with q independent of p to be suitably specified later, reduces equation 
(1) into: 





W = 0 


If we specify q * 2x in equation (4) we obtain: 


(4) 




W = 0 


(5) 


5 


i.e., the Colomb wave equation with fractional angular momentum j 
whereas if we specify q *» 4iT in equation (4) we obtain: 


w 55 + 


_ 1 . il . 4 

4 + § + 


-(« 


H 


( 6 ) 


i,e., Whittaker's equation with independent solutions the M or W func- 
tions with parameters it, ± 3/4 and argument 4/ !,Ti 2 . It thus follows 
from equation (6) and (3) that the two independent solutions f and g of 
equation (1) are: 

( ) - <4ix) 2 ' n^M iT ^ m (4iT!l 2 )...., m - | (?) 

where M are the Whittaker M-functions and T and n are defined by equa- 
tions (2) . Using the properties of the Whittaker functions it is eas- 
ily shown that for t and n real the functions f and g are also real 
with f being even and g odd .functions of n» whereas both f and g are 
even functions of T. The solutions in equation (?) were found earlier 
by us (see Reference [7]) and used in some preliminary studies on wave 
propagation through a linearly sheared flow. 

Series forms for f and g may be easily obtained from the well-known 
series expressions of the Whittaker M-functions. The following is one 
such form: 



n = 0,2,4, . . . for f 
n - 1,3,5, . . . for g 


( 8 ) 


where a^ is determined by the 3-term recurrence relation: 
n 


(4t) 


_ i ; / v 

n = n(n - 3) ( a n-2 ’ a n-4 j , n > 3 


(9) 
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with 


a^ « 1 a, ■ 0 

o i 

o 

a^ * - 8T a^ * 1 

The above series forms were obtained by expanding in series the ex- 
ponential part of the Whittaker M-functions and are in agreement with 
the series obtained earlier by both Kuchemann [2] and Graham and 
Graham [4] although the expressions given by these authors are un- 
necessarily cumbersome. But the expressions in equation (7) for the 
solutions of equation (1) enable one to obtain even batter series-type 
expressions for f and g that may not only be faster converging but also 
more suitable in certain applications. These may be obtained from equa- 
tion (7) using the (series) expressions of the M-functions without ex- 
panding in series the exponential part j nd using the fact that M. (iz) 

it § 

is a real function for real T , m and z. One thus after some algegra ob- 
tains the following forms: 



i- + 2m 
n 2 


cos 


(2tn 2 ) ^ A n (4-nf 


n 


n=0 


+ sin(2Tn 2 ) B n ( 4T, T 2 ) 

n=0 


where the coefficients A and B„ are given by: 

n n 


'a * 

n 

SC 

T - ( 2 " ■— + m) 


A •, 
n-1 

„ B n 


+ (2~~ — 4 m) T 


B„ i 
l n-l J 


( 10 ) 


( 11 ) 
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with m ■ ? j and with the upper sign for f and tha lower for g, and 
*o " 1,1 ®o " °* 

Tha asymptotic expansions for f and g (as well as those of the de- 
rivatives of f and g) may be obtained using Olver's method (6). It 
should be stressed that one is interested in asymptotic expansions as 
T + 00 holding uniformly in n when n ranges over unbounded regions as 
well as asymptotic expansions holding for unbounded Twri+^orO, 
i.e., expansions describing the asymptotic behavior of f and g as func- 
tions of both T and n* The point is not one of mathematical rigour 
only, but also of necessity since the cases n *► 0 or * correspond to in- 
compressible flow or high Mach number limits respectively and t + 0 or 
00 correspond to the long or short wavelength limits > moreover, the cases 
(Tn) - finite characterize the vortex sheet and the incompressible layer 
limits as we will see later. Such uniform asymptotic expansions of the 
Whittaker M- functions have been derived by Skovgaard (9] following 
Olver's method [8]. The results suitable to the present study are in 
terms of the modified Bessel functions of the first and second kind X 2m 
and * 2 m +i* Few of ^he vaulting uniform leading term expressions for the 
asymptotic expansions for f and g for T and r) real we will need later are: 

(i) Leading asymptotic expansion terms for n ** 00 » T-finite: 


-|±mjT 1 

^2-^m ^-r-n T 


[ cos [ - T in (4TT) 2 ) 


Vg/ * m + iT )l » 

+ 2 m 2 + i - J (IT 2 m)j j 

with m « j and 5 * arg J T ( ~ ? m + ix) J 


(12) 
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(ii) Leading asymptotic expansion terms for T -*■ » or ^ + 0. 


0' n (i7TTi? 


fsin h 


cos h 


<4TX) " 4ty <4TX) 

cos h A sin h 


B <x) cos h 
+ « .in h (4TX) 


(13) 


where 


and 


4 2 


|n|/ n 2 - i - An (| n | + /n 2 - ”) 


-i 

32 




M 


5 n 




3 3 

(n 2 - i ) 2 


> (14) 


From the above expressions we may also obtain the following limit 

forms for f and g that we will also need later: 

2 

As n -*■ 00 f T ■* 0 but Tn - finite we obtain from equation (12) : 


^ n 


cos j^2tn 2 - J (l T 2m)] 


(15) 


As ri ■> 0 , x 00 but In - finite we obtain from equations (13) and 
(14): (Note that Aim X = “H and Aim B (x) = 0 or (2x) -1 depending on 


v n+o 

whether n < 1 or r| > !•) 


n-H) 


fv sin h cos h 

K (4th) (4tn> - (4xn) 

g> cos h sin h 


(16) 


* •**'»•*■*»* » •-'•■'MS ’ 


r:yrt3BSXX - % * 


As T '► «■», n > o we obtain from equation (13) : 



sin h 


cob h 


(4tx) 


( 17 ) 


T)>a above asymptotic expressions differ from some corresponding 
expansions obtained by Kuchemann (2) and Pridmore-Brown (31. This 
was expected as mentioned in the introduction since Kuchemann [2] es- 
sentially seeking a high Mach number (or large q) expansion neglected 1 
2 

compared to q in the last term of equation (1) which is tantamount to 
letting T * 0 in Whittaker's equation (6), with the result that the solu- 
tions to equation (6) then are the special Whittaker functions 
2 

M o+ 3/4 (4iTV) ) which are proportional to the Bessel functions 
2 

^±3/4 ) as was obtained by Kuchemann [2) by solving di- 

rectly the approximated equation (1) as mentioned above. But this 
clearly implies that the high Mach number, or large n# solution in terms 
^±3/4^“2tn ) is only valid for T -► 0, i.e. , for low frequencies w. 

In actuality the correct high Mach number asymptotic expansion to lowest 
order is given by equation (12) for all frequencies whereas the low 
frequency/high Mach number limit is given by equation (15) and is ob- 
tained from equation (12) as t *> 0, q + 00 but xq 2 -finite. it is more- 
over easily seen that equation (15) constitutes essentially the leading 
term of the asymptotic expansion of the Bessel functions J ( (-2tq 2 ) 
for large q. The reason for Kuchemann 's [2] incorrect expression 
lies of course in the fact that it is not sufficient to have q 2 »l for 
neglecting 1 in the last term of equation (1) but one should compare 
the term to be neglected to all the other terms of equation (1) which, 
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after some algebra results into the additional necessary condition, 
a -► 0 , .1 . e . , low frequencies or large wavelengths, 'l'ho loading term of 
our asymptotic expansion for T » differs also from that given by 
Fridmore-Brown [3] . It should be pointed out that* Pridmore-Brown [3) 
did not obtain a uniform asymptotic expansion but by applying 
hanger's method [10] he obtained only a nonuniform asymptotic 
(leading) term which is not valid for g -► 0. Our expansions are uni- 
form with respect to both a and n valid in the whole g-plane with the 
exception of an arbitrary neighborhood around the points g * 1. These 
points correspond to points at which the wave normals are parallel to 
the mean flow there. For the propagation problem, these points define 
the location of caustics where rays forming ray tubes converge on a 
line (or a surface) , For the stability problem, as pointed out by 
Betehov and Criminals [6], at these points the pressure fluctuation 
equation changes from elliptic to hyperbolic but the points g 2 » 1 are 
not essential singular points of equation (1) . The behavior of the so- 
lution around these points may bo easily obtained following Skovgaard 
(9) and Olvcr [B] who have obtained uniform asymptotic expansions 
in terms of Airy functions. The leading terms given in equation (13) , 
(16), and (17) are, in particular, uniformly valid in the region g -► 0, 
which, for the propagation problem corresponds to the incompressible 
case and for the stability problem to the solution in the so-called 
"critical layer". 

3. Limiting Cases 

In this section we will examine three limiting casus of our solu- 
tions to equation (1) which correspond to well-known problems of plane 
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wave propagation in a stratified wind and to stability of vortex sheets 
which exemplify the significance of the parameter T and the variable n» 
The limiting cases depending on the values of T and r) are shown in 
Table 1 below: 


TABLE 1 


l 

n *► o 

(a *► <») 

ln| > InJ > 0 

| n | *► «> 
(M * w) 

T ■> 0 

w 

~ *► 0 
b 

w 

a *► < » , £ 0 

but a •> 0 
b 

Incompressible 
Vortex Sheet 
(Kelvin-Helmholtz) 

p + 0 | 

Compressible 
Vortex Sheet 
(Miles, Ribner) 

M 00 and £ *► 0 

(xp) -finite 
Low frequency-lligh 
Mach number 
Finite Shear Layer 

T > > 0 

Critical Layer 

Finite Compressible 
Shear Layer 
(Present Case) 

M ** ® and = finite 
b 

T *► oo 

— CO 

b 

or a 

a *► » , ^ 0 

but (tr)) -finite 
Incompressible finite 
Shear Layer (Rayleigh) 

0) . m 
— oo 

b 

Geometrical Acoustics 
Limit of a Compressible 
Finite Shear Layer 

M •+ oo and —■ *► °» 
b 

1 


We will discuss first the case a 00 which corresponds to the in- 
compressible shear layer and then we will consider the limiting cases 
0 end ~ *+■ 08 , i.e. , the low and high frequency limits which charac- 
terize the compressible vortex sheet and the geometrical acoustics of 
the compressible finite layer respectively. 

3.1. Stability of the Incompressible Shear Layer 
For the finite shear layer, letting a ■+ 00 results in n 0 and 
T 00 # for finite —, and moreover 4'rr) = — - az is finite. In such a 
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case, since our asymptotic expansion is uniform as ^ 0 the solutions 

of the stability equation are those given by equation (16) . Indeed it 
may be easily verified that the two independent solutions of the stabil- 
ity equation for inviscid incompressible fluctuations in a linearly 
sheared flow, i.e. , of the equation 


p nn - t p n ' < 4T|2p * 0 

are exactly those given by our asymptotic forms in equation (16) . 

Moreover, the uniformity of our expansions may also be exhibited 
from the limiting forms of the equation for the eigenvalues. For the 
geometry shown in Figure 1, it can be shown that the eigenvalue equa- 
tion is: 


9 ( n i) 


0 *K) 


- n x f n (n 1 ) 

9 n( n i) 0 

0 f nK) 

■) n (n 2 ) 4 t A - n 2 


(19) 


where f and g are the two independent solutions of the pressure pertur- 
bation equation (1) and T and n have their previous meaning. (See sec- 
tion following equation (2)). 

Using in equation (19) the expressions for f and g given by equa- 
tion (16) , we obtain, after some straight forward though lengthy alge- 
bra, the following expression for the phase speed c. 
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c 1 
b “ 7 


1 

•i 4 |k -■»)* - 1 ci -•,) (* - «- 2a(, ^’)] J 


( 20 ) 


Hie above expression , for z^ ■ 1, z 2 » -1 and b * 1 yields the 
better known special eigenvalue form (Betchov and Criminale 16]). 

c 2 « 1 -i+-i~(i, e “ 4a ) (21) 

“ 4a ' ' 


The limiting form of the incompressible vortex sheet may be easily 
obtained from the above equation (20) in the limit b *► « and z^ and 
z 2 *► 0 but b 2 ^ *► U^, bz 0 *► U 2 , the two constant velocities respectively 
on the opposite sides of the vortex sheet; the result is 


c 


U 1 ♦ U 2 


+ i 



( 22 ) 


which for U * -U “1 reduces to the well-known result c » ± i charac- 
terizing the Kelvin-Helmholtz instability of the incompressible vortex 
sheet. Actually the same result above may be directly obtained not as 
a limit of equation (20) but directly from the general eigenvalue equa- 
tion (19) in the limit n ■+> 0, t ■+• 0 which again exemplifies the uni- 
formity of our asymptotics. In such a case it follows from either equa- 
tion (8) or equation (16)* that f, g and their derivatives f^ , g^ 
acquire the forms; 


Note that the power series equation (8), and the asymptotic forms, 
equation (16), of f and g are identical to second order, i.e., as 
q 0 the uniform asymptotic expansion equation (16) has a second order 
contact with the power series equation (8) . 
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( 23 ) 


* * i - \ (4xn) 2 + o(n 4 ) , « n ^-(4T) 2 n + o(n 3 ) 

9 * j (4in) 3 + —■ (4xn> 5 + o(n 7 ), g n * (4x)V + o(n 4 ) 


Substituting the above limiting expressions in the eigenvalue equation 
(19) and keeping lowest order terras by going to the limit n + 0, x -*■ 0 
we obtain directly the eigenvalues: 



(24) 


which is precisely the result in equation (22) derived previously from 
equation (20). 

3.2. Stability of the Compressible Vortex Sheet 
(Long Wavelength Approximation) 

This case obtained in the limit x + 0 (or ^ •* 0) with T1 finite. 
The series expressions given by equation (8) are to lowest order in Xs 


f^l - <4x) 


„ . . (4t) 

g n + 


' + IL.) -f q/t4\ c «• „ ,_4 


(x**) , f n ^ -(4x) tf (n + n ) + o<t ) 


2 (£-£)♦ 


n/ _4. „ „ ,2 (4x) 

0(x ) , g^ ^ 3ri + — — 


(.*-*)♦ 


0(x 4 ) 


(23a) 


Using these values in the eigenvalue equation (19) we obtain: 



( 25 ) 


The above equation (25) may be easily reduced to a 4th order equation 

for c to yield the eigenvalue expression obtained by Miles [11] and 
others: 
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-%**:** * ' 


'*" 1 " " "" f* *'* '^T: 9^^' >- 5 f^ y h, mu 


I ■ ^i V " 2 ' * ? |m 2 * 4 - 4 /m 2 + 1 [ 2 (26) 

U 1 ‘ U 2 

with M * ~. 

a 

3.3. Plane Wave Propagation Through A Compressible Shear Layer 
(Short Wavelength Approximation) 

This case is usually handled via the so-called Geometrical Acous- 
tics approximation where for the pressure perturbation p(r,t), one as- 
sumes a solution of the form 

nit k n I* . V i (ut-k0) 

p(r,t) * P Q (r,t)e 

where the amplitude p Q (r, t) is assumed to be a slowly varying function 
(of r snd t) and by letting k be large. Finally, by expressing all the am- 
plitudes of the perturbations in series of (ik) -1 one finds that to the 
lowest order in (ik)” 1 , 0 is determined from the so-called eikonal equa- 
tion whereas the amplitude p q may be determined from the energy equation 
as follows (Blokhintsev [12]): 


V0 | 2 - (1 - $ • V0) 2 


(28) 


and 



(29) 


In both equations above the undisturbed and local sound speeds have 

been assumed to be the same and in the energy equation (29) in addition 

to neglecting constant multiplicative factors the partial derivative 

with respect to the time t has also been neglected compared with the di- 
* 

vergence term. For the case of a plane wave incident from a homogenous 


* 

This is justified either in the case of a stationary process, such as 
for a time harmonic field, or on the basis of the original assumption 
of quasistatic perturbation amplitudes. 


half-space (z < 0) at an angle 4> with the z-axis upon a linearly sheared, 

homogeneous medium (z > 0) with velocity U « bzH(z)e where H is the 

x 

Heaviside function, the exact solution of the eikonal equation for the 
phase function 0 is known to be (Kornhauser [13]). 

e « sin 4) x ± ~ [n / X] 2 - 1 - An (n + /n 2 - 1 ) (30) 

The amplitude P Q (r) of the pressure perturbation is not usually 
given in the standard literature but may be obtained from the reduced 
energy equation (29) by observing that the quantity 1 - ft • V0 io, from 
the eikonal equation (28) equal to ± | V0 | » ±n . 

It follows then that 


P Q (r) 'v. 


(n 2 - i)' 


The above results in equation (30) and (31) obtained from the two separate 
equations (28) and (29) of Geometrical Acoustics actually may be obtained 
simultaneously from our uniform asymptotic expansion in equation (13) . 
Since the x-dependence of the pressure perturbation p(r) is 
we may write from equation (27) and (13) in the limit k •* °° or T -> °°: 

9 * Aim -&&- 
k '*’ 0 p o (r) 


ik sin <}) + Jin Jlim /f 

r-K» lg 


J Aim sin h 

l-Kn COS h 


= x sin + —■ An Aim sin h . . 

T-K» cos h 

x sin (j) + |^ r|*4 2 - 1 - An (n + 
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or since 4 t * ~ sin $ * sin <}> in this case, the result in the Geo- 
metrical Acoustics limit equation (30) follows. Moreover, the ampli- 

2 -W 

tudo of the pressure fluctuation r|(q - l) is obtained immediately by 
inspection of our asymptotic expansion, equation (13) . In fact the 
leading term of the asymptotic expansion in equation (13) yields the 
additional information that all the terms of the formal short wavelength 

2 «L 

expansion have the common term n(n - 1) , a result not available from 

the formal Geometrical Acoustics theory (Kornhauser (13) ). 

Finally it should be pointed out that the Geometrical Acoustics 
limit in the example of this section has been obtained from the uniform 
asymptotics of the solutions of equation (1) in the limit t •+• *>, i.e. , 

*► 00 and not in the formal limit xn 2 ^ 00 , according to a criterion that 
has been suggested in the literature (Felsen and Marcuvitz (14)). 
Clearly the later criterion cannot be true for regions arbitrarily close 
to the origin r) = 0 whereas our asymptotic solution is valid through the 
point r) - 0. This observation indicates that care must be exercised in 
the application of the above mentioned (nonuniform) criterion that has 
recently been applied to the solution of certain Aeroacoustics problems 
in the limit of Geometrical Acoustics (Candell (15)). 

2 

4. Physical Meaning of T)» T and 4Tr) 

2 

The physical meaning of f), T and 4xr) which characterize the solu- 
tions f and g of the pressure perturbation equation (1) follow from the 
definitions in equations (2) and the meaning of K which depends on the 
problem at hand. Defining by and the disturbance and relative 
mean flow vector Mach numbers and letting be the unit vector in the 
direction of the parallel mean flow we may write 
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plan* wav* propagation 
stability 

rectangular duct mode* 

( 32 ) 


Tbua r| ii a relative Hach number measure, i.e., it is the parallel to 
the mean flow component of the disturbance (vector) Mach number Mj, 
relative to the relative mean flow Mach number M f . In addition, for 
propagation of a plane wave at incident angle - j < $ < + y from a 
homogeneous half-space <z < 0) on a parallel flow half-space <z > 0) 
it is easily shown that 


n 


l 

sin 4> 


( 33 ) 


where <$> ■ 4>(r) is the (local) angle that the wave normal makes with the 
z axis. 

£i) 

The quantity g acquires the meaning of a characteristic disturbance 
Strouhal number, i.e., 


S - (Shear Layer Thickness) x Disturbance Frequency 
Relative Mean Flow Speed 


Zh> n ID 

bz b 


( 34 ) 


Ihe quantity t is a measure of the disturbance Strouhal number with 
respect to the disturbance Mach number, i.e.. 
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A w „ co/b s Di sturbance Strouhal Huwbgc 

4T " b K * I/K " *► . *► * Component of the Disturbance Mach 

M d *£ Number parallel to the relative 

Mean Flow 


( 35 ) 


2 

Finally the argument 4xn of the Whittaker M-functions also has in 
the case of plane wave propagation the elegant meaning of the (local) 
disturbance wavelength with respect to the relative refractive index 
change. This may be shown as follows: For propagation of a plane wave 
incident from a homogeneous half-space (z < 0) at an angle “ \ S < + 

with respect to the z-axis, on a homogeneous medium (z > 0) with speed 

^ ^ ^ 

U * U(z)e x , it is easy to show that the wave normal unit vectors e n are 

independent of x, i.e, that all wave normals of a given z-stratum are 

parallel. Thus, defining an index of refraction n by 

n * . — i * - — » — - (36 

1 + . t 1 + M 8in * 

x n 


where is the local angle of the wave normal and the z-axis and using 
equation (33) and Snell's Law it follows that 

n ** 1 - M sin 4 > q *= r) sin <J> o * nr) ain <(> (37) 

and one then may write: 

I? n 

Local Disturbance Wavelength o 

Relative Refractive Index Change * |y n | 

n 

sin <p 0 

k o sin $ _ ak o n _ w n 2 m 

~ ^ ^ sin $ " b sin * b sin = 

n a T o 

* F sin n2 s 4Tn2 <38) 
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toJ=l 


This interpretation of the argument of the Whittaker functions of 
the solutions f and g of equation (1) also implies that the necessary 
condition tor the applicability of Geometrical Acoustics suggested in 
the literature (Felsen and Marcuvitz [14]) and recently applied to aero- 
acoustics (Candell (IS)) is not only not a sufficient one, but 
also unnecessary and unnecessarily restrictive. This condition is usu- 
ally given as (Felsen and Marcuvitz (14); Candell (1S))> 



« 1 


(39) 


which, in view of equation (38), implies that 4xn 2 » 1. But in effect 
in Section 3 (3.3) we obtained both the exact solution of the eikonal 
equation and the correct amplitude of the pressure perturbation from our 
uniform asymptotic expansion in equation (13) in the limit x *► 09 indepen- 
dently of the value of r) (i.e., of z or upper fluid Mach number). In- 
deed from equation (13) we see that by comparing the first two terms of 
the asymptotic expansion we easily conclude that the second term is 
negligible compared to the first provided that 4tJx| » 1 which yields 
the criterion equation (39), i.e., 4tn 2 » 1* only in the additional 
limit |n | + 00 . But also since our expansion in equation (13) is uni- 

V 

formly valid for 4 + o we conclude the following two conditions for our 
case (nonzero incidence angle 0 and finite velocity profile slope b ) t 

4th 2 » 1 or ~ « 1 (40) 

z 4 

Y + 0 or \z + 0 (41) 
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Condition (40) above has been interpreted and used ae a relaxed 
criterion for the applicability of Geometrical Acoustics with the impli- 
cation that far enough from the x-axis (large a) the characteristic 
wavelength X need not be unduly small. This interpretation is false in 
view of the second condition in equation (41) which implies that the 
further one is in the far field the shorter the wavelength should be to 
satisfy that condition, -diefia observations indicate that care should 
be exercised in applying the necessary criterion in equation (39) of 
Geometrical Acoustics without examining in detail the behaviour of the 
solutions of the perturbation equations for the particular velocity pro- 
file of the problom at hand. 

5. Applications 

In this section we indicate* how the general solutions of equation 
(1) may be effectively used for certain problems arising in plane wave 
propagation and in the stability of a finite thickness shear layer with 
a linear velocity profile. 

5.1. Plane Wave Propagation 

One of the problems relating to plane wave propagation through a 
finite thickness shear layer with a linear velocity profile is that of 
the existence of "resonances" and the existence of what, in analogy to 
optical wave propagation, may be called Brewster angles. Resonances 
and Brewster angles correspond to infinite and zero values of the 


Here only the general outline of the methodology and certain relevant 
results are given. The details and complete results will be presented 
in another place. 
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reflection coefficient respectively. It is known for instance (Ribner [16], 
Miles [11]) that the vortex sheet, which may be thought as the exact 
limit of long wavelength wave propagation, exhibits both resonances and 
Brewster angles , In contrast short wavelength wave propagation through 
a finite thickness shear layer, i.e., the geometrical acoustics limit, 
does not exhibit any resonance but only a continuum of Brewster angles. 

This is shown in Figure 2. The problem of resonances was touched 
briefly by Graham and Graham [4] , who were only able to show that a 
sufficiently thin, but nonzero thickness, shear layer has no reso- 
nances; they were though unable to draw any conclusions for a finite 
thickness shear layer. The reason for the latter was that Graham and 
Graham essentially obtained and used the series representation, equation 
(8) , of our solution of the pressure perturbation equation (1) . Using 
our solution in terms of Whittaker functions given by equation (7) it 
may be shown that the reflection coefficient R (for the amplitude of the 
pressure perturbation) is given by: 

R* - (CtD ) 2 t42) 

(A ?B) 2 v (C f D) 2 


where 

A - 7F [y°>V u ' V 1,9 n 10) ] 

B * 4T /n 2 - l /nj - 1 [f (l)g(O) - f (0) g (1)3 
C = /n 2 - 1 [f(l)g n (0) - f n (0)g(l)] 

D = / n 2 - 1 [f n (l)g(0) - f(0)g n (l)] 


( 43 ) 
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The upper signs in equation (42) hold for Hj, > 1 and the lower 

signs for < ~ 1» in both cases |r)| >1. In equations (43) we used the 

notation 0 and 1 in the arguments of f and g and their derivatives with the 

understanding that 0 designates evaluation at n ■ H “hi „ * ■ and 

o 1 z>*0 sin 0 

1 designates evaluation at n “ H, * nl * ■— - v - M , i.e., at the 

1 1 z=z^ sin 0 

two edges of the shear layer. (See Figure 1 with z 2 ■ 0.) 

Equation (42) is valid for - y < 0 < + ^ with the upper signs 
holding for the regime of ordinary reflection (r) > 1 and R < 1) and the 
lower signs for the regime of the so-called amplified reflection 

f\ a 

(r^ < - 1# R >1). (For the total reflection regime { R | * 1, 

7T TT 

-l<r,^< + l , - j < 0 < + j .) The conditions for the existence of 
resonances and Brewster angles become: 


A + B»C + D*0 


(44) 

with A - 3 ft 0 

or C ■ DM 

or both. 

Using the Wronskian expression W(n) for our solutions f and g, it 
may be shown that 


AB + CD * -W(0)W(1) 





1 


(45) 


It follows from equation (45) that in general A / B and C / D. 

In order to assess whether the conditions for resonance in equation 
(44) are possible we consider the thin shear layer case, i.e., T small 
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but not zero. Using our expressions in equation (23») we find to lowest 
order in t: 


A + B * 4lM 


" Vii> + A* - l /nf - i {% + n 2 + n c n 1 )J ( 46 ) 


To lowest order in T the right-hand side should be evaluated at the vor- 
tex sheet values of n Q and for which to lowest order in x the second 
condition for the existence of resonances, i.e., C + D a 0 is fulfilled. 
Using the known properties of the vortex sheet solution, i.e.. 


n 2 + n? ~ n 2 rh and n r). » 1 - J m 2 ~+ 1 

o 1 o '1 o 1 


we obtain from (46) : 


A + £ = - (4t) (2M 3 ) 


(47) 


It follows from (47) that A -f B = 0 for a nonzero Mach number only if 
T = 0, i.e., only for the vortex sheet case. The conclusion is that 
the conditions for the existence of resonances equation (44) are incom- 
patible for small but finite T, i.e., for the sufficiently thin shear 
layer and that layer, in contrast to the compressible vortex sheet ex- 
hibits no resonances and no Brewster angles. This is evident also from 
Figure 2 where the reflection coefficient has been evaluated numerically 
as a function of Mach number for various values of T and for a fixed in- 
cident angle of 30° for- which the corresponding compressible vortex 
sheet has two resonances and one Brewster angle. It is also seen in Fig- 
ure 2 that for finite nonzero x the resonances and Brewster angles of 
the corresponding vortex sheet disappear even for very small values of X. 
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5.2. Stability to Long Wavelength Disturbances 
It is known that the incompressible vortex sheet is unstable for all 
wave numbers a whereas the introduction of either finiteness (finite in- 
compressible shear layer) or compressibility (compressible vortex sheet) 
lead to ranges of a or M for which the motion is stable. This is shown in 
Figure 4. The question that naturally arises is then whether the intro- 
duction of both finiteness and compressibility will yield increased or 
any ranges of ci for which the motion is stable. For small a, i.e., long 
wavelengths, this question may be readily answered by the use of our 
series solutions for f and g of the pressure perturbation equation given 
by equations (23a) . For small wave numbers ct or small T using these ex- 
pressions to lowest order in T the eigenvalue equation (19) for the nor- 
malized disturbance phase speed c/a may be shown to be 

x 4 - 2Mx 3 + (m 2 - 2) x 2 + 2Mx - (e + M 2 ) =0 (48) 


where x = ~ , M is the upper fluid Mach number (see Figure 1 with z_ = 0) 

3L * 

and e is the small parameter 


e 



% + 




, M 5* 2 /T 


and 


e is to be evaluated at the vortex sheet value of (~j as determined in 

v 

section 3,2. The allowable solutions of the above equation are given 


exactly by: 


x = 


_ 1 
2 


[m ± /] 


M 2 + 


4-4 


n 


(49) 
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And either by examining the solutions of the 4th order equation (48) by 
use of Rouche's Theorem (Copson £17 1 ) or the algebraic expression for 
the solution, equation (49) above, one easily concludes that for all non- 
zero Mach numbers M (excluding M * 2 J~2) there is always a value of — 
with positive imaginary part. A slight modification of the eigenvalue 

equation yields for M 2 a 5th order equation for x * — and again 

£1 

one may show that c acquires always a complex value with positive imag- 
inary part. 

Thus the long wavelength properties of the shear layer are drasti- 
cally different from those of the compressible vortex sheet and in fact 
contrary to the compressible vortex sheet the thin shear layer is un- 
stable to all Mach numbers for small wave number disturbances. 

6. Conclusions 

In this study we have considered the characterization of inviscid 
fluctuations in a compressible linearly sheared, but otherwise homoge- 
neous parallel two-dimensional flow. The behavior of the cross-flow 
part of the fluctuations were found to be governed by essentially 
Whittaker's equation with the pressure fluctuations being characterized 
by the functions. 

n 1/2 » iT , T 2 (litn 2 ) . 

2 

where M is the Whittaker M-f unctions and r) r T and 4tr) admit the fol- 
lowing simple interpretations 
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ri ■ Disturbance Mach number component in the mean flow direction 
relative to the mean flow Mach number 

T m Disturbance Strouhal number 

Disturbance Mach number component in the mean £ low direction 

4xn 2 ■ Local Disturbance Wavelength 

Relative Refractive Index Change 

The known solutions to a number of other parallel flow problems may 
be obtained as limiting cases of our exact solutions. Such cases in- 
clude the compressible vortex sheet (T ■+ 0) , the incompressible vortex 
sheet (t *• 0, n 0 ) , the incompressible linearly sheared layer (t + 

H ■* 0, (xr|) - finite) and the short wavelength or Geometrical Acoustics 
limit of the compressible linearly sheared layer (t 00 or ^ -*> 0) . 

Finally the exact solutions we have obtained for equation (1) en- 
able us to study exactly the compressible shear layer and to answer 
some as yet unanswered questions such as those pertaining to the exis- 
tence of resonances and Brewster angles and to the stability of such a 
layer. For the shear layer with a linear velocity profile resonances 
and Brewster angles do not exist except in the limits t ■+• 0 (vortex 
sheet) and x ■* °° (geometrical acoustics) . Moreover for small wave num- 
ber disturbances the linear shear layer is unstable for all Mach numbers M 
Thus contrary to the compressible vortex sheet, which is stable for 
M > 2 / 2, one may not regard the compressible vortex sheet as an adequate 
model of a thin shear layer for all relative Mach numbers of the uniform 
flows bounding the shear layer. 
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LIST OP CAPTIONS 


Figure 1; 
Figure 2: 


Figure 3: 


Figure 4: 


The linear velocity profile shear layer geometry. 

o 

Reflected energy flux/incident energy flux (- R ) versus 
h or upper fluid Mach number M for T » 0 (vortex sheet) 
and x - 00 (geometrical acoustics) . Shear layer geometry 
as in Figure 1 with z 2 ■ o. Incident acoustic plane wave 
at an angle 6 ■ + 30° (clockwise) with the positive z-axis. 
Reflected energy flux/incident energy flux (* R ) versus n 
or upper fluid Mach number M for various values of T. 

Shear layer geometry as in Figure 1 with ■ 0. Incident 
acoustic plane w ?»ve at an angle 0 ■ + 30° (clockwise) with 
the positive z-axis. 

Regions of stability of vortex sheets and shear layers with 
a linear velocity profile. 
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